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A general constructive procedure is presented for analyzing magnetic instabilities in two- 
dimensional materials, in terms of [predominantly] double nesting, and applied to Hartree-Fock 
HF+RPA and Gutzwiller approximation GA+RPA calculations of the Hubbard model. Applied 
to the cuprates, it is found that competing magnetic interactions are present only for hole doping, 
between half filling and the Van Hove singularity. While HF+RPA instabilities are present at all 
dopings (for sufficiently large Hubbard U), in a Gutzwiller approximation they are restricted to a 
doping range close to the range of relevance for the physical cuprates. The same model would hold 
for charge instabilities, except that the interaction is more likely to be q-dependent. 
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I. INTRODUCTION 

The charge and magnetic instabilities of one dimen- 
sional materials are well understood in terms of Fermi 
surface (FS) nesting. Here it is shown that for two- 
dimensional materials, features in maps of the bare sus- 
ceptibility can be understood in terms of a series of FS 
'nesting curves', and the dominant instabilities are gen- 
erally related to double nesting features. The analysis 
shows how to locate these nesting features in momentum 
(q) space, often providing analytical expressions. As an 
application, the full evolution with doping of the leading 
magnetic instabilities for several families of cuprates is 
presented, both in the conventional Hartree-Fock (HF) 
plus RPA (HF+RPA) and in a Gutzwiller approxima- 
tion (GA+RPA) calculation. The analysis provides a 
pseudogap candidate and makes the surprising predic- 
tion that the 'checkerboards' seen in scanning tunneling 
microscopy (STM) studies are not the same phase as the 
'stripes' in La2- :r Sr. E Cu04 (LSCO). 

The present results should find extensive utilization. 
First, for weakly correlated two-dimensional systems the 
HF+RPA results provide an essentially complete solution 
to the nesting problem, as long as the interaction U is q- 
independent. The GA+RPA extends the results into the 
intermediate coupling regime. In both cases, the mag- 
netic instabilities are determined by zeroes of the Stoner 
denominator, 

l-tf e //Xo(q,w = Q). (1) 

Here for cuprates U e ff is the Hubbard U in HF+RPA, 
and a more complicated q-dependent object Uga(<i) m 
the GA+RPA calculation. Thus, the leading HF+RPA 
instability is simply associated with the maximum of 
the bare susceptibility xom — max q XoCo.; 0), while the 
leading Gutzwiller instability can be shifted by the q- 
dependence of Uga(<1)- It will be clear that the same 
analysis can be extended to any two-dimensional mate- 
rial. 



The cuprates appear to be in an intermediate coupling 
regime where the Gutzwiller results can be expected to 
provide a good approximation to the phase diagram at 
T = 0. Thus, in the cuprates, recent quantum Monte 
Carlo (QMC) 1 and 'quasiparticle-GW (QP-GW)P calcu- 
lations have been able to reproduce experimental ARPES 
spectra of optimally and overdoped cuprates, starting es- 
sentially from LDA bands and calculating the self-energy 
self consistently. In the QP-GW approach, the self en- 
ergy is calculated as a convolution over a renormalized 
RPA susceptibility. Not only is the low-energy disper- 
sion rcnormalization reproduced, but also the 'water- 
fall' effectP, which represents the dressing of low-energy 
quasiparticles by (mainly) spin fluctuations. Extension 
of these results to the hole-underdoped regime will re- 
quire identification of the phase or phases responsible for 
the pseudogap. Since these are most likely to be incom- 
mensurate density wave or 'stripe' phases, the QMC cal- 
culations have a severe problem of limited q-resolution, 
while the QP-GW calculations are ideally suited to han- 
dle this. Such an approach has had considerable success 
with electron-doped cuprates, where the leading instabil- 
ity is commensurate (vr, tt). (We set the lattice constant 
a = 1.) 

The paper is organized around the Stoner criterion 
as follows. Sections 2-3 describe the calculations of the 
zero frequency bare susceptibility Xo(q)i an d determin- 
ing its maximum Xom- In Section 2 we introduce the 
concept of nesting curve, associated with the nesting cri- 
terion q = 2kp, where is the (anisotropic) Fermi 
wave vector, and we demonstrate that xom is gener- 
ally associated with a double nesting. We apply this 
concept to the cuprates, taking the hopping parameters 
from tight-binding one-band fits to the LDA dispersions, 
and show that over most of the doping range there is 
a unique Xom, with a surprising electron-hole symme- 
try. In Section 3 we extend this analysis to the more 
complicated hole doping regime for doping x between 
half filled and the doping of the Van Hove singularity 
(VHS). Here competing magnetic phases are found, and 
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the dominant phase is sensitive to material parameters, 
being different for different cuprates. Section 4 describes 
the GA+RPA technique and introduces the correspond- 
ing U e ff — Uga(cl)- Section 5 presents the resulting 
Gutzwiller magnetic phase diagrams. It is found that the 
phase diagram of La 2 _ a; Sr. E Cu04 (LSCO) is distinct from 
that of most other cuprates, but that for all cuprates, 
using a bare Hubbard U = 8t, the paramagnetic state is 
unstable at the GA level for all relevant dopings, includ- 
ing the overdoped regime. A discussion is presented in 
Section 6, and conclusions in Section 7. 



II. SUSCEPTIBILITY PLATEAUS 

The present calculations are based on tight-binding 
parametrizations of typical dispersions for single-layer 
cuprates, including models of La 2 _ x Sr ;r Cu04 (LSCO), 
Nd 2 -xCe x Cu0 4 (NCCO), and Bi 2 Sr 2 Cui0 6 (B12201), 
Table I in Appendix [Aj For NCCO, Bi-2201(2), and 
LSCO(2), the parameters arc based on a tight-binding fit 
to the LDA bandP, while sets Bi-2201(1) and LSCO(l) 
are direct fits to the experimental bands. In all cases, 
k z dispersion is neglected, approximating the cuprates 
as two-dimensional. 

The real part of the susceptibility Xo(l> w = 0) is dom- 
inated by a series of plateaus, with the largest suscepti- 
bility systematically shifting from one plateau to another 
as a function of doping. Figure [T] shows that these and 
related features dominate x'o over the full doping range, 
and further reveals a striking quasi-electron-hole symme- 
try of the evolution. However, instead of being symmet- 
ric about x = 0, the evolution is symmetric about the 
doping of the Van Hove singularity (VHS), xvhs, where 
the Fermi energy coincides with the VHS, Ep = Eyas, 
and the density of states (DOS) diverges logarithmically. 
Figure [TJa) shows x'o calculated along the high symmetry 
axes of LSCO, for the full electron-doping range n = 1—x 
from 1 to 2, with an extension to the hole doping of the 
VHS (blue curves). Near half filling the susceptibility is 
dominated by the well-known platearP near [tt, tt) (red 
and blue curves), which peaks at the VHS. As electron- 
doping increases, the peak shifts to a second plateau near 
(tt, 0) (green curves), then to a third near T for a nearly 
full band (brown curves). The same sequence is followed 
for hole doping, Fig.[IJb). 

In addition to these plateaus, there is an additional F- 
centered feature, which is prominent near the doping of 
the VHS, but has largely disappeared by half filling. This 
feature will be referred to as the antinodal nesting (ANN) 
plateau. Bi2201 displays the same plateaus, with the 
same quasi-electron-hole symmetry, Fig. []Jc,d). A key 
difference is that the (tt, 7r)-plateau of Bi2201 is convex 
over a considerably wider doping range than in LSCO, 
before turning concave near the VHS. This has a conse- 
quence that the ANN peak dominates over a wide doping 
range in Bi2201, but only in the immediate vicinity of the 
VHS in LSCO. At the VHS, the bare susceptibility di- 
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FIG. 1: (Color online.) (a) Susceptibility x'o( u = °) for LSCO 
as a function of q for a series of dopings from xvhs = 0.207 
(blue dashed line) to x — —0.99 (brown dashed line). [Blue 
curves are for hole doping, x > 0, others for electron doping, 
x < 0.] Note the evolution from a plateau near (it, tt) (blue 
and red lines) to one near (tt, 0) (green) to one near V (brown). 
Band parameters appropriate to LSCO(2) (Table I), (b) Same 
as (a), except for hole dopings x = 0.207 (red dashed line) 
to x — 0.99 (brown dashed line). Note similar evolution of 
plateaus. (c,d) Similar to (a,b), except for Bi-2201 (2 in Table 
I). Note overall similarity, except for curvature near (tt, 7r). 
Dashed curves correspond to x = (c): -0.99 (brown) or 0.43 
(blue); (d): 0.44 (red) or 0.99 (brown). Curves are generally 
spaced by Ax = 0.05, except for (1) higher density near points 
of rapid change (e.g., the VHS), (2) Ax = 0.1 near top and 
bottom of band, and (3) end points at x — ±0.99. 

verges logarithmicly both at T (where it is equal to the 
DOS) and at (7r,7r), but the latter divergence is quite 
weak and not apparent in the numerical calculations of 
Fig. 1. Despite the apparent complexity of these evolv- 
ing susceptibility patterns, the plateau evolution can be 
understood in detail, with analytic formulas for the po- 
sitions of all dominant peaks. 

The generic evolution of the plateaus from (tt, tt) — > 
(tt, 0) — > (0, 0) can be understood with reference to the 
'nesting curves', Fig. [2} For the generic case of two Fermi 
surface segments, a nesting curve can be defined as the 
locus of all points q = kpi — kp2, where kpi is a point 
on the ith FS, FSi, with the restriction that when FSi 
is shifted by q it is tangent to FS 2 . For the cuprates, 
there is usually just a single FS section [an exception is 
given below in Section III, Fig. 6], and the nesting curves 
simplify to plots of 

q = 2k F , (2) 

for any Fermi momentum k F , Fig. [2j Since the point 
(tt, tt) lies on horizontal and vertical planes of reflec- 
tion symmetry for the susceptibility, the nesting curves 
must be supplemented by their reflections about the lines 
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FIG. 2: (Color online. )(a) Hole-like Fermi surface (solid) 
shifted by the double nesting vector q = (0.38, 1.62)7r/a 
and folded back into the first BZ (dashed), (b) Electron- 
like Fermi surface (solid) shifted by the double nesting vector 
q = (0.58, l)n/a and folded back into the first BZ (dashed). 



q x = 7T, q y = 7T. Fig. [3] shows three sets of nesting curves 
corresponding to three hole dopings for LSCO (dispersion 
1), x = 0.41 (red curves), 0.62 (blue), and 0.79 (green). 
Frames (b-d) show the corresponding susceptibility maps 
Xo(q), and it can be seen that the ridges in x'o are ex ~ 
actly given by the nesting curves. [For convenience, the 
x = 0.41 data is replotted in Fig. 3(a). Note that the 
susceptibility maps are plotted over the range of q x , q y 
between and 7r, whereas the nesting curves are plotted 
over the wider range to 2ir.] Furthermore, the dominant 
peaks in x'o correspond to the intersection of two nest- 
ing curves. By drawing the original and q-shifted Fermi 
surfaces (FSs), Figs, ga), (b), it can be seen that this 
overlap corresponds to the simultaneous nesting of two 
different sections of FS. Hence the term 'double nesting'. 

As shown in Figs. [2^c) , (d) this kind of 'double nest- 
ing' can originate from either the scattering between 
points on opposite segments (cf. example in Fig. [2^c)) 
or adjacent segments (cf. example in Fig. [2](d)) of the 
FS. In the example shown in Fig. [2^c) we denote the 
scattered states on the FS as Pi = (5±,n — 6\\) and 
P2 = (— 7T + <5|| , — 6±) which yield the points on the nest- 
ing curves qi = 2Pi = (25_l,27t-25||) = (2<5_l,-25||) and 
q 2 = 2P 2 = (-2tt + 2£||,-2£l) = (25\\,-26±). Thus in 
this case 'double nesting' (qi = q 2 ) occurs when 5± = S\\ 
which generally can only be fulfilled for points near the 
antinodes of hole- like FS's. These are thus referred to 
as antinodal nesting (ANN) features. More common is 
the situation sketched in Fig. [2ld) where we consider the 
scattered states Pi = (6±,w—6\\) andP 2 = (S±, — ir+6\\). 
In this case the 'double nesting' condition qi = q 2 can 
only be fulfilled for 6\\ = (which is trivial since initial 
and final states of the two scattering processes are iden- 
tical) or <J|| — 7r/2. The latter condition implies that this 
kind of 'double nesting' generally affects the nodal states 
and leads to scattering vectors q close to (ir/a, it J a), gen- 
erating the (7T, 7r)-plateau. 

A slightly different point of view might help clarify the 
role of double nesting. Tangency of two Fermi surface 
segments ['nesting'] means a stability of FS overlap, in 
that the nesting persists if one displaces the nesting vec- 
tor in a particular direction. Then double nesting means 
tangency along two surface segments, so nesting persists 
if one displaces the nesting vector in two (rather than 
one) directions. 

The origin of the plateau transitions is now apparent: 
as hole doping increases, the FS shrinks to a small pocket 
near T before disappearing. The nesting curves shrink 
in a similar fashion, but with a doubled radius, since 
q = 2kp. Thus the dominant overlap shifts from near 
(7r,7r) at the VHS (Fig. 2(b)) towards (tt,0) (Fig. 2(c)), 
and finally towards T, in Fig. 2(d), thereby explaining 
the plateau evolution. For electron doping, the Fermi 
surface ultimately closes at (tt,tt), leading to the same 
sequence of plateau transitions, as illustrated in Fig. |1J 
where the nesting maps are directly superposed on the 
susceptibility curves. In the doping range relevant to the 
cuprates, the physics of NCCO is dominated by the (%, ir) 




FIG. 3: (Color online.) (a) Nesting maps q = 2k F for LSCO, 
dispersion 1, and x = 0.41 (red curves), 0.62 (blue), and 0.79 
(green). Corresponding susceptibility Xo( w = 0) as a function 
of q for a series of hole dopings from Ef [x] = -0.16 [0.41] (b), 
-0.3 [0.62]] (c), and -0.5 eV [0.79] (d). In all figures, whites 
are largest x s > blues are smallest. 

plateau, Fig. [4^, and red curves in Fig. 1(a), which shrinks 
to a point at the end of the 'hot-spot' regime. 

In all cases, the dominant peak lies along a high sym- 
metry axis, and the doping dependence of its position 
can readily be found from the dispersion 

E(k) = -2t[c x (a) + Cy{a)] - Wc x (a)c y (a) 
-2t"[c x (2a) + c y (2a)} 
-4t'"[c x (2a)cy{a) + c y {2a)c x {a)] , (3) 

where 

Ci(aa) = cos(aki<i) , (4) 

and a is an integer (or half-integer). For this dispersion, 
the VHS is generally at (tt, 0), or E VHS = 4f - 2t" . The 
positions of the peaks in Fig. (3jh,c) and Fig. |4^b,c) lie 
along the zone boundary (q x a = tt) with 

q y a — 2 arccos[±-\A 2 — c — b], (5) 

with b = (t- 2t"')/[At"] and c = [E F - At") /[At"]. For 
electron [hole] doping the peak is exactly at (tt, 0) when 
E F = +[-]2(f - 2t"'). Beyond this point, the peak is 
at q y — 0, q x given by Eq. pj] with b — (t — [+]2t' + 
2t"')/[A(t"-[+]2t"')] ande = \E F -[+]2(t-2t"')]/[A(t" - 
[+]2t"')]. Since a peak at T corresponds to ferromag- 
netism, the above reproduces the common finding that 
a nearly empty or nearly full band tends to be ferro- 
magnetic. The ANN peak qi = q2 in Fig. [2^l satisfies 
cos(q y a) = [E F - 2t']/[2{t' - 2t")}. Thus the peak sus- 
ceptibility is generally associated with double nesting. 
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FIG. 4: (Color online.) Susceptibility Xo( w = 0) as a function 
of q with superposed nesting maps q — 2k F for NCCO, for a 
series of electron dopings (x < 0) from x [Ef] = -0.26 [0.10] 
(a), -0.37 [0.20] (b), -0.52 [0.35] (c), and -0.65 [0.50 eV] (d). 

The only exception we have found to this is a tendency 
to remain commensurate for a finite doping range about 
high symmetry points such as (7r,7r) or T — (0,0)J^- a 
form of Van Hove nesting-^. 

III. COMPETING PHASES AND 
FERROMAGNETISM 

While the above sequence is completely generic, hold- 
ing for all cuprates and being electron-hole symmetric, 
the additional features associated with the VHS are con- 
siderably more variable. Fig. [5] illustrates the low-hole 
doping regime for three dispersions proposed for the 
cuprates. At half filling, all three have a peak suscepti- 
bility on the (tt, tt) plateau, associated with conventional 
staggered antiferromagnetism (SAF). At low doping the 
peak is either commensurate or (jr, tt — 5) incommensu- 
rate, with 8 increasing with doping. For finite hole doping 
a new feature emerges, a roughly +-shaped peak along 
the zone diagonal, at (tt — S, tt — 6) . From the nesting 
curves, it can be seen that the peak is associated with 
nesting of the flat sections of the bands near (tt, 0) - hence 
the name antinodal nesting (ANN) - but the largest sus- 
ceptibility lies along the zone diagonal, where both (tt, 0) 
and (0, tt) nesting occur simultaneously, Fig. [2p. This 
ANN feature has an interesting relation with 'hot spot' 
physiciP. A 'hot spot' is a point of the Fermi surface 
which is simultaneously on the antiferromagnetic (AF) 
zone boundary [diagonal of the paramagnetic Brillouin 
zone]. From Fig. [2jc) it becomes apparent (cf. large 
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FIG. 5: (Color online.) Susceptibility maps for Bi2201 (a-d), 
LSCO(l) (e-h), and LSCO(2) (i-1) at a series of increasing 
dopings. x [E F \ = 0.12 [-0.20] (a), 0.20 [-0.25] (b), 0.31 [- 
0.30] (c), 0.40~ xvhs [0.33] (d); 0.19 [-0.12] (e), 0.23 [-0.13] 
(f), 0.29 [-0.14] (g), (x VHS ~ 0.33), 0.37 [-0.15] (h); 0.0 [-0.05] 
(i), 0.14 [-0.10] (j), 0.20~ xvhs [-0.111] (k), and 0.21 [-0.1112 
eV] (1). 

dot) that the image of the AF zone boundary in the ex- 
tended q-map intersects the 'nesting curves' exactly at 
the points of 'double nesting'. But the image of the AF 
zone boundary gets folded in the q-map onto the zone 
diagonal, T — > (tt,tt), thereby generating the usual ANN 
feature. Thus all the diagonal ANN peaks in q arise from 
hot spots in k. As the doping of the VHS is approached 
the hot spots move toward (7r, 0) and the nesting curve 
moves to T. Since the susceptibility at T is equal to 
the DOS, it diverges at the VHS, thereby controlling the 
magnetic instability. 

The appearence of the ANN peak leads to a com- 
petition between two different kinds of magnetic order, 
and the doping evolution of the susceptibility maps di- 
verges. For most dispersions studied, including the left 
and central columns of Fig. [5j the ANN intensity grows 
and the dominant peak changes discontinuously from the 
SAF plateau to ANN nesting. For the Bi-2201 disper- 
sion (left column) this ANN peak becomes dominant at 
about x—0.2, and evolves smoothly to T at x=0A. If 
one unfolds the nesting curve, one sees that this happens 
exactly at the doping of the VHS, when the FS passes 
through (7T, 0). The central column, corresponding to an 
extreme dispersion proposed for LSCO to enhance antin- 
odal nesting, displays a much more complicated nesting 




FIG. 6: (Color online.) (a) Susceptibility map for LSCO(l) 
(h) with nesting curves superposed. Red curves: q = 21cf 
nesting for both T-centered barrel and (tt,0), (0, 7r)-centered 
pockets. Blue lines = barrel-pocket nesting; green line = 
inter-pocket nesting, (b) Corresponding Fermi surface map 
showing four q-shifted pockets, labelled as in (a). 



map over a limited doping range close to the VHS, but 
also displays a dominant peak at T exactly at the VHS. 
This map will be explained below. Finally, for the right 
column, corresponding to a more conventional dispersion 
for LSCO, the ANN peak is weaker than the SAF peak 
except in the immediate vicinity of the VHS. In this case, 
the dominant susceptibility peak remains commensurate 
at (tt, tt) until the doping of the VHS, jumps to T at the 
VHS, then jumps back and smoothly evolves to (tt, tt — 5) 
incommensurate [this sequence can also be followed in 
Fig. 1(a), (b)]. Note further, in Fig. 5(i,j), that there is a 
wide doping range where the susceptibility peak remains 
commensurate, at (tt, tt). [Since (tt, 7r) lies along several 
mirror planes, the corresponding susceptibility is in gen- 
eral a local maximum or minimum.] 

The complicated nesting map of the middle column 
near the VHS is explained in Fig. [6] The dispersion is 
such as to produce an 'extended VHS'P which first in- 
tersects the r — > (7r,0) axis at a point (tt — 6,0). For 
larger hole doping, the FS has two sheets, one a squar- 
ish barrel centered at T, the other a pocket centered at 
(tt, 0), Fig. [6^b) . The intrasheet q = 2kp nesting maps 
are shown by the red curves in Fig. |6ja). In addition 
to the one associated with the barrel FS, there are two 
overlapping segments associated with the pockets, but 
translated by the q folding to the vicinity of T. As in the 
other frames of Figs. 2-4, these nesting curves exactly 
match some of the ridges seen in the susceptibility. How- 
ever, there are additional ridges, associated with nesting 
between two different FS segments. Thus, the blue curves 
in Fig. |6ja) represent barrel-pocket nesting - that is, the 
locus of the q-vectors needed to shift the pocket until it 
is tangent to some point on the barrel FS. For instance, 
the points labelled 1, 2, 3 represent translations equiva- 
lent to those shown in Fig.[6]jb). The new nesting curves 
can easily be found numerically, by requiring that the two 
FS sections have a common tangent at the point of oscu- 
lation. Similarly, the green curve near (tt, tt) represents 



inter-pocket nesting, with point 4 illustrating one partic- 
ular nesting vector. It can be seen that the full collection 
of nesting curves explains all of the ridge-like features 
seen in the susceptibility map, and in particular allow 
the determination of the points of maximal susceptibil- 
ity, except for the above-noted commcnsurability effects. 
Note however, that despite these complications, the sus- 
ceptibility peak moves to T in a finite doping range about 
the point where the antinodal electron pocket shrinks to 
zero width. 

The above discussion can be summarized: in the dop- 
ing between half filling and the VHS, a new susceptibil- 
ity peak arises, associated with nesting of the antinodal 
parts of the large FS. Two kinds of behavior are found: 
when the ANN is dominant, the peak evolves to T at the 
VHS. On the other hand, for LSCO(2) with small t', the 
ANN peak is inherently weaker than the (jr, tt) peak, in 
which case the susceptibility remains commensurate at 
(it, 7r) from half filling almost to the VHS, then smoothly 
develops a (ir, ir — S) incommensurability. Yet even here, 
in the immediate vicinity of the VHS, the ANN peak 
at r becomes dominant in a very limited doping range. 
Hence, for most dispersions the susceptibility will have 
strong FM fluctuations near the VHS, while for other 
dispersions the fluctuations remain mostly AFM. 
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IV. GUTZWILLER CALCULATION AND U, 



In the GA+RPA calculatio n 1 10 ! 11 !, the electronic para- 
magnetic ground state energy is calculated in the G A 13 1 14 1 
and then expanded to second order in the on-site and in- 
tersite magnetic fluctuations in the spirit of Vollhardt's 
Fermi liquid approach^. Response functions are com- 
puted using linear response in the presence of small ex- 
ternal field™ One obtains RPA like susceptibilities but 
with strong vertex corrections. As shown in Appendix [B] 
longitudinal and transverse susceptibilities are trivially 
related and lead to the same Stoner criteria as required 
by spin rotational invariance. In terms of a tensor bare 
transverse susceptibility (cf . Appendix [B] and Appendix 
C) 
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with Ek q = ek +q + e^, the GA+RPA dressed suscepti- 
bility Xq is given by 
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(7) 



Here the ratio of dressed to bare dispersion is given by 
£ k/ e k = Z i with the Gutzwiller renormalization factor 
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FIG. 7: (Color online.) Susceptibility maps (a) x?ii 
and (c) X22> along with derived quantities (d) Ei and (e 
as a function of q for Bi2201 with x = 0.15. 



(b) X°i2, 

E2, 



which depends on the GA double occupancy variational 
parameter D and the doping x. The interaction matrix 
isP 



v+- 

q 



VM Q J' 



(9) 



(8) 



which is defined in Appendix |B| 

While Eq. [7] is a tensor equation, it can be expanded 
into the form of Eq. ^ with xo — z oXii and U e ff = Uga, 

U ga - (^Vq + 2M q £ 1 + Ml{Et - El) X0 /zl)/zl (10) 

with Ei = X12/X11, El = xVXii- The z factors ap- 
pear now in Uga since we want to use the bare xo in 
Eq. [l] All correlation effects are incorporated in the def- 
inition of Uga which therefore can be viewed as a vertex 
corrected interaction term in the magnetic p — h chan- 
nel. Note that Uga is not strictly a pure interaction 
term, but is weighted by kinetic terms which enhance its 
q-dependence. 

Figure [7] compares x?i ; X12 > X22 > an d the derived E\ , 
E2 - the results are fairly insensitive to doping or band 
parameters. The sharp structures associated with nest- 
ing show up only in xli, while the other %'s are smooth, 
and the E's contain only a weak structure imposed from 
Xii- Hence Uga remains a smooth, weakly varying func- 
tion of q, and the instabilities are controlled by the peaks 
in xo- Figure [8^ a) shows how Uga varies as the bare U 
is increased. At small U, Uga —> U and the conventional 
HF+RPA is recovered. Again the results do not depend 
strongly on doping (Fig. 8(b)) or dispersion. The general 
trend is that at large U Uga tends to saturate, and the 
large-q components of Uga are reduced most strongly, 
tending to favor instabilities nearer T — (0,0), sugges- 
tive of ferromagnetic domains. Fig.|8|should be compared 
with Fig. 2 of Vilk and Tremblajff^I who find a similar 
saturation. Since these authors employ a very different 
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FIG. 8: (Color online.) (a) Variation of Uga with U for 
x — 0.01. (b) Variation of Uga with x for U = 2Ubr- Uga 
varies with q, and the figure shows q = T (red solid line), 
(n, 0) (blue dashed line), (n/2,n/2) (brown dotted line), and 
(•7T, 7r) (green dot-dashed line). The q-points along the zone 
diagonal, (ir, 0) and (n/2, n/2) have nearly identical behavior. 

perscription for including vertex corrections, the similar- 
ity of our results gives considerable further confidence 
to the trends we find. In addition since the Gutzwiller 
approximation has a simple interpretation in terms of ki- 
netic energy suppression due to correlation our results 
shed further light on the physical meaning of the vertex 
corrections. 



V. GUTZWILLER MAGNETIC PHASE 
DIAGRAMS OF THE CUPRATES 

Figure [9] shows our main result, the magnetic phase di- 
agrams of LSCO and Bi2201 calculated in the Gutzwiller 
GA+RPA approach, based on the Stoner criterion, Eq. 1. 
The phase diagram for Bi2201 is fairly generic, and we 
have found similar results for NCCO, Bi2212 [neglecting 
bilayer splitting] and SCOC, while the LSCO phase dia- 
gram is limited to the dispersion parameters of LSCO(l). 
Thus, for electron doping the phase diagram is very sim- 
ple, dominated by simple AF order with q very close to 
(■7T, 7r), while for hole doping there is a competition be- 
tween the SAF and ANN order, except in LSCO, where 
SAF order dominates. It is possible that ANN order is 
relevant to the 'checkerboard' phase seen in STM studies 
of several cuprates, as discussed briefly below. 

Results for HF+RPA are similar, but since U e fj hi 
Eq. 1 is not bounded {U e ff = U), there would be mag- 
netic instabilities for all dopings, with transitions to the 
(tt,0) and T plateaus, as in Fig. [I] In contrast, in the 
GA+RPA calculation Uqa saturates, and magnetic or- 
der exists only in a limited doping range in the cuprates, 
primarily in the range where the (tt, tt) plateau exists. 
However, this doping range encompasses the full dop- 
ing range of relevance to cuprate physics - including the 
overdoped regime - if a doping independent U ~ 6 — 8i 
is assumed. Thus, to have a magnetic quantum criti- 
cal point near optimal doping, one must go beyond the 




X 



FIG. 9: (Color online.) Gutzwiller magnetic phase diagrams 
for (a) LSCO(2) and (b) Bi2201(2). Dashed lines indicate 
metastable states - extensions of the condition UautzX = 1 f° r 
one phase beyond the point where another phase has become 
unstable. 

one-band GA+RPA approach. The limited doping range 
for the magnetic phases has also been seen in a d = oo 
studypH Analogously to HF+RPA, GA+RPA determines 
the transition line at which the mean field state becomes 
unstable, but the energies are evaluated at mean field 
level. Fluctuation contributions to the energy will most 
likely lower the energy of the paramagnetic (disordered) 
phase further reducing the stability range of the magnetic 
phases. 



VI. DISCUSSION 

There have by now been countless calculations con- 
cerning nanoscale phase separation or stripe phases in 
cuprates, and this is becoming a paradigm in many other 
correlated electronic systems as well. We believe, how- 
ever, that the present results are unique in providing 
a systematic phase diagram, covering the full doping 
dependence^ and all possible q- vectors for realistic band 
dispersions, and indeed finding that the dominant insta- 
bilities may be different in different cuprates. 



A. Extension to Charge Instabilities 

While the above has provided a thorough analysis of 
possible magnetic instabilities, nothing has been said 
about competing instabilities in the charge channel - 



charge density waves (CDWs) or 'stripes'. Of course, 
one possibility is that charge order arises as a secondary 
effect following the spin order, with 5p oc (SM) 2 . The 
case where the charge instability is primary will be some- 
what harder to analyze. That is because the corre- 
sponding Stoner factor analysis is likely to require an 
extension beyond the Hubbard model, and hence in- 
volve greater uncertainty in the choice of t/ e //(q). For 
charge instabilities, the corresponding Stoner denomina- 
tor is nothing but the (zero frequency) dielectric constant 
e/e = l+V(q)xo{ c l, w = 0), where e is a background di- 
electric constant and V(q) is a Coulomb potential. When 
V(q) is taken as either the Hubbard U or a long-range 
Coulomb interaction,^ e is always positive and there is 
no instability. 

Moreover, from the theory of dielectric stability, it is 
known that e cannot fall in the range between and 1 
[equivalently eT x < 1], and an instability e = must 
be approached through negative values of eP^ While 
a purely electronic instability could still arise via in- 
clusion of local field effects, the most natural situation 
arises when electron-phonon coupling is included. A sim- 
ple s wave instability will be suppressed by the on-site 
Coulomb repulsion as shown in Ref . [TD] but one can still 
have instabilities in other channels. To the extent that 
the matrix elements which determine the channel have a 
smooth behavior one will have a Stoner denominator sim- 
ilar to Eq. (ll with Uj^f ~ 9 1 : A) / 'e e ^t 2 phl where g is the 
electron-phonon coupling parameter, Dq the correspond- 
ing bare phonon propagator, a bare phonon fre- 
quency, and e e an electron-electron dielectric constant.^ 

Thus, since the instability is also controlled by a Stoner 
factor, the results of the present paper will also apply for 
charge instabilities. This is consistent with the common 
expectation that CDWs and spin density waves (SDWs) 
are controlled by the same nesting instabilities. Stated 
differently, the Stoner criterion is a formal expression of 
the idea that P/K > 1 for an instability, where P, the in- 
teraction energy, is represented by U e f / while the kinetic 
energy K is measured via Xo ■ The kinetic energy is 
closely tied to the band structure, and carries important 
material-specific information. In contrast, the potential 
energy is fairly featureless, leading to smoothly varying 
U e ff(q). Thus, the locations of the instabilities are con- 
trolled by peaks in xo, and a q-dependent U can only 
shift the dominant instability between two peaks of com- 
parable height. In the remainder of this section, when we 
discuss comparison to experiment, it will be seen that the 
present model has correctly determined the dominant q 
values, but that in several cases experiment points more 
toward CDWs than SDWs. 



B. Stripes vs Checkerboards 

There have been a number of recent hints that the 
'stripe' order in underdoped La^^A^CuO^^, A = Sr 
(LSCO) or Ba (LBCO) is not the same phenomenon as 



the 'checkerboard' order found in Sr 2 Cu0 2 Cl 2 (SCOC) 
and Bi 2 Sr 2 CaCu 2 8 (Bi-2212). Thus, resonant soft 
x-ray scattering experiments find evidence of charge 
order in LBCO (stripes)^ but not in Ca 2 Cu0 2 Cl 2 
(checkerboards)^, while evidence for time reversal 
symmetry breaking has been found in YBa 2 Cu307_a 
(YBCO), but not in LSCO^S. The present results sug- 
gest a connection between the ANN and checkerboard 
phases, in that the periodicity of the latter also scales 
with the antinodal nesting vectoi^, and a further con- 
nection between the (-7T, n) plateau and stripes. Thus it is 
quite interesting to observe that the ANN phase is virtu- 
ally absent for the LSCO (2) dispersion, where evidence 
for conventional stripes is strongest. 



C. Conventional and VHS nesting 

In ID systems the susceptibility diverges at q = 2kp 
at all dopings, due to perfect (flat band) nesting. The 
present results generalize this to 2D systems: the domi- 
nant instability is generally associated with double nest- 
ing, but since xo m general does not diverge, a finite 
coupling U is required to drive an instability. An ex- 
ception is VHS nesting^, for which xo has a logarithmic 
divergence. 

At the VHS, there are two competing instabilities, AF 
at (7r, 7r) and FM at T, corresponding to the magnetic 
branch of the 50(8) phase diagram of the cuprates, and 
these are responsible for the commensurate pinning near 
(ir, 7r) and T, respectively.^ Doped away from the VHS, 
these instabilities evolve into the two dominant peaks 
of the susceptibility, and can be considered as 'general- 
ized VHS nesting'. As such, they dominate the mag- 
netic physics of the cuprates over the full doping range, 
coextensive with the limits of the (tt, 7r)-plateau in the 
susceptibility. 



D. VHS and FM 

Recently, Storey et al^ proposed that the generic be- 
havior of high-T c cuprates could be understood if the 
pairing interaction (or pairing energy cutoff) falls off 
rapidly near the VHS. A recent experiment!^ does indeed 
confirm that x op t scales with T c , but x op t/xvHS < 1- 
The present results suggest that most cuprates will have 
strong FM fluctuations near the VHS, which are in- 
compatible with simple d-wave superconductivity. Em- 
pirical evidence for FM fluctuations has been noted 
previously^. While some previous calculations^ have 
found evidence of ferromagnetism near a VHS, others^ 
have suggested that FM instabilities were unlikely in 
competition with incommensurate susceptibility. The 
present calculations confirm that a dominant FM sus- 
ceptibility should be present near the VHS. 
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E. Limitations of Present Approach 

1. k z - dispersion 

The advantage of two-dimensional materials is that it 
is straightforward to display the susceptibility maps and 
nesting curves. In 3D, the curves become nesting surfaces 
and the susceptibility maps are 4-dimcnsional. One can 
speculate that the dominant susceptibility peaks corre- 
spond to triple-nesting points. For quasi-2D materials, it 
should be possible to analyze a series of cuts perpendic- 
ular to the (weakly-dispersing) z-axis. 



TABLE I: I. Band Parameter Sets 
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2. Nanoscale phase separation 

In principle, the present results could provide informa- 
tion on nanoscale phase separation (NPS) as well. One 
model of NPS is that there are different instabilities as- 
sociated with particular dopings [e.g., half filling and the 
VHS], and that these two phases are more stable than 
uniform phases at intermediate dopings. In this case, 
the two end phases could still be described by nesting 
maps, only at particular dopings. For instance, in the 
HF+RPA analysis of LSCO(2), the FM phase is stable 
only very close to the VHS, and could lead to NPS with 
a second phase at the undoped insulator. 

3. Away from the Instability Threshold - Towards Strong 
Correlations 

The Stoner criterion determines which q value is most 
unstable, and the minimum U needed to drive that in- 
stability. However, as U increases above threshold, the 
q of the ordered phase may shift. Thus, at half filling, 
when a full gap can be formed, q will no longer be de- 
termined by best nesting, but by the largest gap. This 
tends to favor more commensurate q values, leading to a 
pinning of q at these commensurate values over a wide 
range of parameters. We find that as U increases, there 
is a first order transition from an incommensurate phase 
with Fermi surface pockets to a commensurate [(tt, it) 
or (tt, 0)] phase that is fully gappedPS Away from half 
filling, Luttingcr's theorem ensures the persistence of a 
Fermi surface, so nesting instabilities should persist over 
a wider range of Vs. 

Furthermore, as we have seen above, the q-value cor- 
responding to the largest /7 e //(q)xo(q) can shift with 
U. Typically, for Uqa the shift is to a smaller q-value, 
associated with an instability towards FM order (in the 
magnetic channel) or NPS (in the charge channel). 

VII. CONCLUSION 

The present results provide a constructive scheme for 
identifying the dominant nesting instabilities for any two 



dimensional material. Clearly, for Fermi surfaces with 
multiple sections, a large number of nesting curves are 
possible, leading to extremely complicated susceptibility 
maps. Nevertheless, the present scheme will automati- 
cally sort out the possible double nesting peaks and follow 
their evolution with doping. This should allow a much 
more detailed understanding of 2D phase diagrams, par- 
ticularly for magnetic phases, where the interaction U 
has negligible q-dependence. 

The good agreement of the GA+RPA calculations with 
more detailed variational calculations at half filling 34 and 
with exact and numerical results in d = oo, 11 and of the 
AF+SC model in electron doped cuprates with experi- 
ments gives us confidence in the model. Accordingly, we 
note the following three points. First, within GA+RPA, 
the paramagnetic state is unstable over the full hole- 
doping range in the cuprates, including overdoped. To 
avoid this conclusion, and restore a Fermi liquid phase in 
the overdoped regime, it may be necessary to include non 
Gaussian fluctuations^ or an additional doping depen- 
dence for U. The origin of any such doping dependence 
lies outside the GA treatment of the Hubbard model. 
Second, the (jr, tt) phase, or its incommensurate exten- 
sion, is unstable against a competing ANN order in [most] 
hole-doped cuprates. Third, there is a material depen- 
dence to the phase diagram, and LSCO may have a very 
different doping dependence than other cuprates. 

The results have been applied to a number of model 
dispersions for the cuprates. A possible pseudogap candi- 
date has been identified and a distinction between stripes 
and checkerboards proposed. These findings will be dis- 
cussed in greater detail in ensuing publications. 
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ern University's Advanced Scientific Computation Cen- 
ter (ASCC). RSM's work has been partially funded by 
the Marie Curie Grant PIIF-GA-2008-220790 SOQCS, 
while GS' work is supported by the Vigoni Program 2007- 
2008 of the Ateneo Italo-Tedesco Deutsch-Italienisches 
Hochschulzentrum . 
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Appendix A: Notes to Table I: Band Parameters 

In Table I, all hopping parameters are given in meV. 
The t"" term for LSCO(l) is the coefficient of a term 
c x (2a)c y {2a) in Eq. [3] NCCO, LSCO(2), and Bi2201(2) 
data sets are fit to LDA band dispersions of the near- 
Fermi level antibonding CuC>2 band, and are appropriate 
for the Gutzwiller analysis. The data sets LSCO(l) and 
Bi2201(l) are taken from fits to experimental ARPES 
data. To the extent that the bands are renormalized by a 
q-independent factor Z, the effective bare susceptibilities 
Z\ can be used in the Gutzwiller analysis, while Xom 
occurs at the same q. 

Recently it has been suggested that the enhanced nest- 
ing, such as observed in LSCO(l), may be associated with 
stripe formation^. 

Appendix B: Gutzwiller plus random-phase 
approximation formalism. 

Here we sketch the GA+RPA formalism and define the 
relevant quantities. For a review see Ref. 1361 

The susceptibility can be computed in the 
longitudinal 10 or in the transverse channel.^ As 



shown bellow both results are equivalent as dictated by 
spin rotational invariance on a paramagnetic (singlet) 
ground state. 

For later use we define the density matrix associated 
with the unprojected Slater determinant, \<p), as 




and as a shorthand we define pij a = pff. 

We start from the spin-rotational invariant Gutzwiller 
energy functional for the Hubbard model defined in terms 
of p1° and the double occupancy in the Gutzwiller vari- 
ational state Di. As derived e.g. in Ref.l37lthe functional 
reads, 

e ga = £ u, (^iU^ \4,) + u A- 

i,j i 

where we have defined the spinor operators 

•'-(2) 

and the z-matrix 



z^ cos 2 ^ + Zii sin 2 ^ ^ [ztf - z^] cos (pi 



§g? fo't - zqj cos ip t z A sin 2 ^ + z iX cos 2 ^ 



tan 2 ^ = ''' "' 



with the z factors given by 

_ yfi Pa + D ^Pu + ^ - A) + y/A(|p» - ^ - Di) 

Zia ~ r~\ 's 1 1 ~s~ z — 



and for clarity spin expectation values are denoted by 

St = P f U> S i = Pu f > S * = (P f u ~ P^)A and p u = 
Pii +Pit- 

In the limit of a vanishing rotation angle <p the z-matrix 
becomes diagonal and the renormalization factors reduce 
to those of the standard Gutzwiller approximation. 

For a homogeneous, paramagnetic system (zj CT = Zo> 
ifi = 0) the expansion of the Gutzwiller energy functional 
up to second order in the particle-hole excitations reads, 



SE GA = 5E Sf + SE- 



ga 



SE, 



GA 
lo 



SE\ 



GA 



(Bl) 



The first contribution is the free fermion part, 



E 

k>kp : k' <kp;(T,CT' 



£k denotes the dispersion of the Gutzwiller quasiparti- 
cles whereas corresponds to the unrenormalized dis- 
persions, i.e. £k = z 2 e k - 



The remaining part in Eq. (Bl| is due to quasiparti- 



cle interactions and separates into contributions from the 
charge (c), the longitudinal (lo) and transverse (tr) spin 
channel. SE GA and SE GA have been derived in detail in 
Ref. EH 

For our purposes we report in the following the expan- 
sion in the longitudinal channel: 
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and in the transverse magnetic channel which has been 
derived in Ref. [37] 
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The relevant fluctuations have been defined as 
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The longitudinal susceptibility describes the Am z = 0, 
singlet to triplet excitations of the paramagnetic state 
while the transverse describes Am z — ±1 spin excita- 
tions. Spin rotational invariance dictates that these ex- 
citations should be degenerate. 

One obtains for these correlation functions 



0,lo 
Aq 



and 
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For the non-interacting Gutzwiller quasiparticles spin- 
rotational invariance is thus guaranteed from the relation 



1 



^q 2 -^-Q 
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(B9) 



<S2? = 



The 11 element of the interaction kernel is given by 
7V q - (zJ) 2 iV lq + 2z z' t 'e (B6) 

with 



(B7) 



k<r 



Here Ubr is the Brinkman-Rice energy. N\ is similar, 
but requires separate averages for different components 
of the energy: 
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in an obvious notation. 

Furtheron we define the susceptibility matrices for the 
bare time-ordered correlation functions both in the lon- 
gitudinal 



This identity is preserved within the GA+RPA when 
we compute the interacting susceptibilities from the RPA 
series 
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since from Eqs. (B3 B5) the interaction kernels are re- 
lated by V zz — 2V tr . Clearly the energies of particle 
hole excitations described by these equations are degen- 
erate. In particular they lead to the same Stoner criteria 
in both channels. 



Appendix C: Derivatives of the Gutzwiller 
approximation z factors 

The expressions involve the following derivatives of Zi, 
in the longitudinal channel 



dz 



dz lo 
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and in the transverse channel 

n,tr (f] = zl( (TS+(t)SZ q (0))o (TS+ (t)f I q (0)) \ 
Xq U N \ (Tf+(t)SZ q (0))o (Tf+(t)f =q (0)) o ) ' 

where a hat has been added to distinguish fluctuations 
(SSg) form operators (<S+). 
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The derivatives in the transverse channel can be related _ 1 ^ 

to those in the longitudinal channel as follows, 2 (1 — 5 — 2D) 2 

i _ dzi^-a i , 28 ( 1 



B 2 z i 



Qp a : "dp a,C! 2 " On the right we have given explicit expressions in terms 

2z , / 1 ^ ^ of 5, the doping measured with respect to half-filling and 



(1 -5 2 ) 2 



1 — 2(5 I -j — 1 ] > the double occupancy D in the paramagnetic state. 
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